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CHAPTER  I 


INTRODUCTION 
PROBLEM  STATEMENT 

This  research  effort  encompassed  an  area  in  need  of  study 
that  had  been  identified  by  the  U.S.  Air  Force  Flight  Dynamics 
Laboratory  (AFFDL)  at  Wright-Patterson  AFB,  Ohio.  The  AFFDL 
Structures  Division  is  responsible  for  obtaining  structural  cost 
estimating  data  which  is  used  in  the  development  of  aircraft  design 
to  cost  estimates.  Cost  data  is  currently  obtained  from  a computer 
program  which  designs  an  aircraft  structure  subject  to  a minimum 
weight  routine.  Since  minimum  weight  design  does  not  inherently 
result  in  an  optimum  cost  design,  the  current  program  routine  may 
produce  excessive  cost  estimates.  In  order  to  upgrade  the  existing 
structural  design  and  analysis  computer  program,  a feasibility  study 
was  needed  to  determine  whether  a series  of  design  algorithms  could 
be  developed  which  would  perform  an  optimum  tradeoff  between  mini- 
mum cost  and  minimum  weight  ( 1 1 ; 13;  19). 
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BACKGROUND 
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During  the  mid  1960s,  the  Department  of  Defense  (DOD) 
utilized  a total  package  procurement  (TPP)  policy  for  major  weapon 
system  acquisitions.  This  procurement  policy  called  for  a single 
point  decision  strategy  (10:4),  Industry  was  awarded  a contract  to 
design,  develop,  and  produce  a weapon  system  for  a fixed  price. 

Under  TPP,  the  cost  parameter  was  often  obfuscated  by  the 
syndrome  of  "nothing  but  the  best. 11  The  TPP  policy  was  eventually 
undermined  by  high  costs,  unprecedented  inflationary  trends,  the 
Vietnam  war,  and  schedule  risks  (5:60).  In  1968,  primarily  as  a 
result  of  excessive  cost  overruns.  Deputy  Secretary  of  Defense 
David  Packard  initiated  a major  study  effort  in  the  area  of  defense 
materiel  acquisition  to  examine  current  policy  and  make  recommenda- 
tions for  improvement  (2).  From  this  study,  the  impetus  emerged 
for  a complete  revision  of  DOD  acquisition  policy.  The  revision 
culminated  in  the  publication  of  DOD  Directive  5000.  1 in  July  1971 
which  formalized  the  procedures  for  acquisition  of  major  defense 
systems  and  established  cost  as  a co-equal  to  design  specifications 
and  system  performance  (23:4).  Specifically  the  directive  required 
that: 

Cost  parameters  shall  be  established  which  consider 
the  cost  of  acquisition  and  ownership;  discrete  cost  ele- 
ments (e.g.  , unit  production  cost,  operating  and  support 
cost)  shall  be  translated  into  "design  to"  requirements. 
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System  development  shall  be  continuously  evaluated  against 
these  requirements  with  the  same  rigor  as  that  applied  to 
technical  requirements  [23:4]. 

The  concept  of  design  to  cost  <DTC)  was  embraced  by  DOD 
in  order  to  use  unit  cost  goals  as  thresholds  for  managers  and  as 
design  parameters  for  engineers  (22:41.  The  DTC  process  was 
formalized  in  1971  by  the  publication  of  DOD  Directive  5000.  28. 

This  directive  established  the  framework  for  DTC  whereby  DTC 
goals  would  prevent  funds  from  being  spent  beyond  the  point  where 
costs  rise  rapidly  for  small  increments  of  increased  performance. 
However,  minimum  performance  requirements  were  not  to  be 
sacrificed  (22:4). 

Experience  to  date  indicates  the  objectives  of  DTC  have  not 
been  met  and  cost  growth  in  DOD  major  weapon  system  acquisitions 
continues  to  be  a problem  (19).  Past  practices  in  the  acquisition 
process  are  continually  being  examined  to  determine  whether  they 
are  still  needed  or  require  modification  (4:1).  A study  of  45  weapon 
system  cost  overruns  conducted  by  the  Comptroller  General  in  1973 
revealed  a cost  growth  of  $19.  1 billion  over  development  estimates 
of  $93.6  billion.  Analysis  indicated  three  major  causes  for  this 
dramatic  growth  of  costs:  improved  weapon  capability,  requirements 
change,  and  inadequate  cost  estimating  ability  (6:23-4).  In  the  U.S. 
Air  Force  A-10A  close-air  support  aircraft  development  program, 
the  DTC  goal  grew  from  $1.  524  million  in  197  3 to  $2,215  million  in 
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June  1976.  Twenty-two  percent  of  this  cost  increase,  $.  151  million, 
was  attributed  to  errors  in  estimating  the  original  goal  (12).  Based 
on  production  of  729  aircraft,  this  would  amount  to  S 1 1 0 million  (1:26). 

If  DTC  is  to  remain  viable,  improved  techniques  for  estimat- 
ing costs  must  be  developed  to  insure  DTC  goals  are  set  realistically 
and  represent  the  best  tradeoff  possible  between  cost  and  perform- 
ance (19).  These  methods  must  be  applied  early  in  the  conceptual 
phase  of  the  weapon  system  acquisition  process.  Analysis  of  cost 
estimating  techniques  performed  on  the  F-16  air  combat  fighter  and 
the  avionics  system  of  the  B-l  bomber  has  indicated  significant  cost 
reductions  can  be  achieved  by  applying  detailed  cost  study  methods 
during  the  early  stages  of  the  conceptual  phase.  In  the  F-16  pro- 
gram, a cost  savings  of  over  $2.  5 billion  was  achieved  by  making 
several  relatively  simple  configuration  changes  (19).  The  Boeing 
Aerospace  Company  cost  studies  increased  the  research,  develop- 
ment, test,  and  evaluation  (RDT&cE)  costs  of  the  B-l  avionics  system 
by  $5.4  million.  However,  there  was  a reduction  from  estimated 
production  costs  of  $105.7  million  (20:13).  These  savings,  in  the 
form  of  cost  avoidance,  were  the  direct  result  of  early  cost  studies 
(19).  Results  of  the  F-16  and  B-l  studies  indicate  that  savings  can 
be  achieved  through  more  accurate  cost  estimating  techniques. 

To  assist  AFFDL  Structures  Division  in  improving  their 
cost  estimating  capability,  the  Vehicle  Design  and  Evaluation 
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Program  (VDEP)  was  implemented  in  1975  (see  Appendix  E).  VDEP 
was  developed  by  General  Dynamics  Convair  Aerospace  Division 
under  contract  to  the  National  Aeronautics  and  Space  Administration 


(NASA)  to  aid  in  the  preliminary  design  of  aerospace  structures.  The 
program  provides  a vehicle  synthesis  capability  that  includes  vehicle 
sizing,  structural  analysis,  and  cost  evaluation  ( 1 6 : 1 ) . 

The  cost  evaluation  portion  of  VDEP  is  driven  by  optimiza- 
tion routines  that  attempt  to  design  a minimum  weight  airframe.  It 
provides  no  automated  search  or  optimization  routines  that  attempt 
to  design  structures  at  minimum  cost.  The  cost  of  a particular  air- 
craft design  can  be  estimated  by  using  the  program,  but  it  is  not 
known  whether  a minimum  cost  design  has  been  obtained  (19). 

A unique  feature  of  VDEP  is  an  interactive  systems  package 
which  allows  dynamic  change  of  program  input,  thereby  enabling  the 
designer  to  perform  tradeoff  studies  using  various  structural  vari- 
ables. While  this  capability  has  potential  for  improving  cost  estima- 
tion, it  is  complicated  by  the  fact  that  there  are  more  than  360,  000 
possible  variable  combinations  (see  Figure  1),  which  can  be  input. 
Since  each  combination  of  variables  costs  $11.35  to  run  through 
VDEP,  this  amounts  to  more  than  $4.  1 million  to  obtain  cost  esti- 
mates with  no  guarantee  that  the  results  will  be  cost  effective  ^9). 

The  complexity  of  the  interactive  systems  package  and  the  high  cost 
implications  of  such  an  approach  emphasize  the  need  for  algorithms 
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to  assist  structural  design  engineers  in  developing  a minimum  cost 
design. 

SCOPE 

The  task  of  accomplishing  a feasibility  study  and  developing 
algorithms  has  been  estimated  to  require  approximately  two  man- 
years,  exclusive  of  reporting  time,  and  over  $100,  000  in  computer 
costs  (19)-  This  research  effort  could  not  accomplish  :he  complete 
task  within  available  time  and  resource  constraints.  Therefore,  it 
was  conducted  to  determine  whether  the  development  of  such  algo- 
rithms was  feasible  for  use  with  VDEP.  The  results  provide  a 
foundation  for  future  research  work  in  the  development  of  algorithms 
to  assist  in  achieving  minimum  cost  designs. 

OBJECTIVE 

The  objective  of  this  research  effort  was  to  determine, 
through  a feasibility  study  and  computer  simulation,  whether  struc- 
tural design  cost  optimization  algorithms  could  be  developed. 

RESEARCH  QUESTIONS 

To  accomplish  the  research  objective,  this  study  addressed 
the  following  questions: 
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1.  Which  variables  or  variable  combinations  should  be  used 


as  primary  inputs  for  developing  algorithms? 

2.  Among  the  selected  variables,  what  relationships  exist 
for  developing  algorithms? 

3.  Can  algorithms  be  developed  for  use  in  VDEP  to  produce 
an  optimum  cost  aircraft  structure  design? 


CHAPTER  II 


METHODOLOGY 
DATA  COLLECTION 


Data  Source 

The  data  necessary  for  conducting  the  feasibility  study  was 
generated  by  the  CDC  6600  computer  system  located  at  the  Aero- 
nautical Systems  Division  (ASD),  Wright-Patterson  AFB,  Ohio.  The 
computer,  together  with  the  VDEP  software  package,  provided  an 
aircraft  design  capability  which  was  used  to  study  the  effects  of 
selected  variables  on  aircraft  cost  and  weight.  ;See  Appendix  B for 
a detailed  discussion  of  VDEP). 

In  conjunction  with  VDEP,  the  DC- 10  data  base  was  selected 
for  the  following  reasons: 

1.  Currency.  The  DC-10  information  represented  the  most 
current  data  package  available  for  use  with  VDEP  (13). 

2.  Accuracy.  Data  could  be  readily  compared  with  an 
operational  vehicle  (13). 

3.  Relevancy.  The  DC-10  is  among  several  aircraft  cur- 
rently being  considered  as  an  Air  Force  Advanced  Tanker-Cargo 


Aircraft  (17). 
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Univers  e 


The  universe  for  data  collection  was  the  DC- 10  aircraft. 
Population 

The  population  was  the  wing  structure  of  the  DC- 10.  iSee 
Appendix  A for  operational  definitions  of  wing  and  structural 
variables . ) 

Samples 

The  samples  consisted  of  a series  of  random  observations 
of  cost  and  weight  based  on  the  average  thickness  of  wing  skin 
(T-bar),  material  type,  structural  type,  load  factor,  and  rib  spacing. 
Figure  2 depicts  the  wing  structural  concepts. 

Variables 

The  variables  selected  for  generating  sample  observations 
were  the  result  of  recommendations  made  by  personnel  of  the  AFFDL 
Structures  Division  (19).  Wing  structural  variables  selected  for  use 
in  the  study  were  integral  blade  (LI),  integral  zee  <L3),  integral  tee 
(L4),  and  separate  jay  (L5).  (See  Figure  2.)  The  material  used  was 
aluminum  (19).  The  structural  variables  along  with  T-bar  (T), 
weight  (W),  rib  spacing  (R),  and  load  factor  iU),  comprised  the 
independent  variables.  The  dependent  variable  was  the  total  cost  of 
the  DC- 10  wing  including  the  costs  associated  with  tooling. 
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manufacturing,  overhead,  and  material  acquisition.  Cost  was  com- 
puted as  a function  of  various  combinations  of  the  independent 
variables . 


Sampling  Plan 

The  sampling  plan  was  developed  in  accordance  with  recom- 
mendations made  by  the  AFFDL.  Their  guidance  provided  the  vari- 
able selection  and  priorities  for  data  gathering  used  in  this  research 
effort.  The  rib  spacing  parameter  was  limited  to  values  between  26 
and  76  inches  with  a cost  of  design  being  produced  for  each  increment 
of  10  inches.  The  T -bar  values  were  generated  by  inputting  load 
factors  in  increments  of  100,  000  from  300,  000  to  600,  000  pounds; 
increments  of  200,  000  from  600,000  to  1,  800,  000  pounds;  and 
increments  of  600,  000  from  1,  800,  000  to  3,600,000  pounds  ( 19). 

Using  aluminum  as  the  material  type,  data  was  collected  as 
the  values  for  load  factor,  rib  spacing,  and  structural  types  were 
varied  within  the  limits  stated  above.  This  generated  a total  of  312 
observations  for  the  resultant  T-bar,  cost,  and  weight  data.  The 
weight  of  the  wing  was  recorded  at  each  observation  for  subsequent 
cost  versus  weight  analysis. 

STATISTICAL  ANALYSIS 

Multiple  Linear  Regression 

Multiple  Linear  Regression  (MLR)  analysis  provides  a 


A 
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highly  useful  analytical  method  for  experimental  situations  where  the 


experimenter  is  able  to  exercise  control  over  the  independent  vari- 
ables (14:215).  Due  to  the  nature  of  the  data  that  was  collected, 

MLR  was  the  appropriate  statistical  model  to  be  applied  to  determine 
whether  there  was  a relationship  between  cost  and  the  independent 
variables,  T-bar,  structural  type,  load  factor,  rib  spacing,  and 
weight.  (See  Appendix  C for  a complete  description  of  this  statistical 
analysis  technique.) 


Partial  Correlation 


This  statistical  technique  was  used  to  determine  whether 
spurious  correlations  existed  among  the  variables. 

Spurious  correlation  is  defined  in  a relationship 
between  two  variables,  A and  B for  example,  in  which  A's 
correlation  with  B is  solely  the  result  of  the  fact  that  A 
varies  along  with  some  other  variable,  C for  example, 
which  is  indeed  the  true  predictor  of  B.  In  this  case, 
when  the  effects  of  C are  controlled,  held  constant,  etc.  , 

B no  longer  varies  with  A [15:303]. 

If  the  fixing  of  one  variable  reduces  the  correlation  between 
two  variables,  then  it  is  inferred  that  their  interdependency  arose  in 
part  through  the  agency  of  the  control  variable  and  if  the  partial  corre- 
lation was  zero  or  very  small,  it  is  inferred  their  interdependence  is 
entirely  attributed  to  the  fixed  variable.  If  the  controlled  partial 
correlation  is  greater  than  the  original  correlation,  then  it  can  be 
inferred  that  the  controlled  variable  is  masking  the  true  relationship 


(7:317). 
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CRITERIA  TEST 


The  results  of  the  statistical  hypothesis  tests  on  the  MLR 
model  were  used  to  indicate  which  individual  variables  provided 
statistically  significant  net  or  marginal  contributions  to  the  explana- 
tory power  of  the  MLR  model.  Only  significant  independent  variables 
were  retained  in  the  model  since  these  variables  exhibited  the  great- 
est impact  on  cost,  the  dependent  variable.  The  retained  variables 
and  their  relationships,  as  specified  by  the  model,  formed  the  basis 
for  determining  whether  algorithms  could  be  developed  for  use  with 
VDEP  to  derive  minimum  cost  aircraft  structural  designs. 

ASSUMPTIONS  AND  LIMITATIONS 
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Limitations 

1.  The  study  was  conducted  using  cost  figures  in  1971 
dollars.  Therefore,  the  conclusions  must  be  viewed  accordingly. 

Any  extrapolation  to  other  year  dollars  without  the  use  of  appropriate 

financial  techniques  will  be  in  error. 

2.  VDEP  was  initialized  to  provide  cost  and  weight  data 

for  the  200th  aircraft  structure.  Generalizations  drawn  to  a struc-  1] 

ture  other  than  the  200th  one  will  not  be  valid. 

I 


CHAPTER  HI 


DATA  ANALYSIS 
INTRODUCTION 

In  this  chapter  the  results  of  MLR  are  analyzed,  as  well  as 
the  results  obtained  from  partial  correlation.  Cost  was  regressed 
wiln  all  independent  variables:  weight  (W),  rib  spacing  (R),  load 
factor  (U),  T-bar  (T),  and  structural  types  (Li).  The  material  type 
used  was  aluminum  (15).  Structural  types  were  integral  blade  (LI), 
integral  zee  (L3),  integral  tee  (L4),  and  separate  jay  (L5).  Integral 
tee  was  encoded  as  a dummy  control  variable  and  therefore  did  not 
appear  in  the  cost  equations. 

Only  the  statistically  significant  independent  variables  were 
retained  from  the  initial  MLR.  Cost  was  then  regressed  with  these 
significant  variables  to  obtain  a new  model  for  further  analysis. 

The  next  portion  of  the  chapter  deals  with  the  MLR  of  cost 
with  the  selected  independent  variables  weight,  T-bar,  and  structural 

types  on  an  individual  basis.  Results  from  these  regressions  are 

o 

analyzed  for  significance.  The  MLR  of  cost  with  all  independent 
variables  except  T-bar  and  again  with  all  but  weight  are  examined  to 
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determine  what,  if  any,  effect  these  two  excluded  variables  have  on 


the  overall  regression. 

Analysis  of  the  results  obtained  from  the  regression  of  cost 
with  all  significant  independent  variables  and  cost  with  selected  vari- 
ables pointed  out  interrelationships  which  necessitated  the  use  of 
partial  correlation.  1 he  res lifts "of  the  partial  correlation  runs  are 
presented  and  explained. 

The  final  section  of  the  chapter  examines  the  results  obtained 
from  plotting  cost  with  weight  and  cost  with  T-bar  for  all  other  inde- 
pendent variable  combinations.  These  graphs  depict  the  effects  of 
weight  and  T-bar  on  cost,  and  the  point  at  which  the  least  cost  is 
found  for  each  variable  combination. 

REGRESSION  ANALYSIS 

Cost  with  All  Independent  Variables 

The  results  of  MLR  for  cost  with  all  independent  variables 
are  shown  in  Table  I.  The  variables  entered  the  regression  equation 
using  the  forward  inclusion  scheme.  The  following  hypothesis  test 
for  the  overall  regression  was  performed, 

H0:RZ  = 0 

H^R2  t 0. 

Computed  F (F0)  was  obtained  using  p-1  |7)  upper  and  n-p 
(304)  lower  degrees  of  freedom  (d.f.  ).  F0  = 3428.  18,  therefore 
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Table  I 

Multiple  Regression,  Cost  with  All 
Independent  Variables 


Variable 

Name 

Bi 

Fo 

W 

Weight 

82.  03 

506. 07  * 

R 

Rib  Spacing 

-181 15. 82 

491.71  * 

LI 

Integral  Blade 

-245121. 52 

21.  39  * 

• U 

Load 

*0.  48 

4. 80  ** 

T 

T - bar 

5371  13. 71 

2.  89 

L5 

Separate  Jay 

-35105. 52 

0.  84 

L3 

Integral  Zee 

-6063. 52 

0.  03 

Cons  tant 

2559270. 88 

- _ 

Significant  at  a = 0.0 1 

5jC  S{C  « 

Significant  at  a = 0.  35 


Cost  = 2559270.  88  + 82.  03  W - 181 15.  82  R - 245121. 52  LI  - 0.48U 
+ 537113.71  T - 35105.85  L5  - 6063.52  L3 

Total  Variation  (TV)  = 125.  12  x 10^ 

Explained  Variation  (EV)  = 123.  55  x 10^ 

Unexplained  Variation  (UV)  = 1.57  x 10^3 

R2  = 0.  9875 

Regression  FQ  = 3428.  18 


P(F>FQ)  was  less  than  0.001.  Since  the  F value  computed  was 
highly  significant,  the  null  hypothesis  (HQ)  was  rejected.  The  con- 
clusion was  that  R2  does  not  equal  zero  statistically.  Therefore,  at 
least  one  B^X^  was  significant,  i.  e. , at  least  one  independent  variable 
had  a relationship  to  cost.  Further  hypothesis  tests  were  conducted 
on  the  regression  coefficients.  The  hypothesis  test  on  the  net  or 


marginal  contribution  of  X = was: 

J 1 

H0:Bi  = 0 

H^Bj  * 0. 

Using  the  conservative  Bonferroni  technique,  an  equivalent  a -level 
was  computed: 

equivalent  a = — - — , where  K = other  independent  variables 
K+  1 

in  the  model.  To  obtain  the  Fc  for  a given  equivalent  a-ievel,  it  was 
necessary  to  use  the  Wonnacott  and  Wonnacott  technique  of  interpola- 
tion (26:600).  Comparing  the  individual  F statistic  to  the  interpo- 
lated Fc;  T-bar  (T),  integral  zee  (L3),  separate  jay  (L5),  and  load 
(U)  were  found  to  be  insignificant  at  a < 0.  01. 

Cost  with  Significant  Variables 

Table  II  shows  the  results  obtained  from  regressing  cost 
with  the  significant  variables.  The  hypothesis  test  for  was  con- 
ducted and  the  null  hypothesis  was  rejected  at  all  levels  of  a > 0.001. 
Similarly,  all  variables'  B^X^  were  found  to  be  significant  for 
a>  0.01,  using  the  equivalent  a technique. 

Cost  with  Selected  Variables 

Table  III  reflects  the  results  of  regressing  cost  with  selected 
independent  variables.  In  both  the  regression  of  cost  with  weight  and 
cost  with  T-bar,  the  models  were  highly  significant.  However,  when 
cost  was  regressed  with  the  structural  types,  L5  \separate  jay)  had 
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Table  II 

Multiple  Regression,  Cost  with  All 
Significant  Independent  Variables 


Variable  Name  FQ 

W Weight  78.72  22648.  06 

R Rib  Spacing  -17870.  38  560.32 

LI  Integral  Blade  -194940.94  37.78 

Constant  --  2515130.  22 

Cost  = 2515130.  22  + 78.72  W - 17870.  38R  - 194940.94  LI 

TV  = 125.  12  x 10i3 
EV  = 123.  52  x 1013 
UV  = 1. 60  x 1013 

R2  = 0.  9872 

Regression  FQ  = 7945.48 

an  insufficient  tolerance-level  and  did  not  enter  the  model.  From 
the  foregoing  analysis,  it  was  concluded  that  structural  type  (Li) 
taken  alone  is  not  a significant  contributor  to  cost. 

Cost  with  All  Variables  Except  T-bar /Weight 

Table  IV  depicts  the  results  of  regressing  cost  with  all 
independent  variables  except  T-bar.  No  unexpected  relationships  were 
found.  The  overall  regression  was  significant  at  a > 0.001. 

Table  V shows  the  results  of  regressing  cost  with  all  inde- 
pendent variables  except  weight.  The  overall  regression  was  signifi- 
cant. Load  (U)  entered  the  equation  first  since  it  had  the  greatest 
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Table  III 

Multiple  Regression,  Cos t with  Selec ted 
Independent  Variables 


Cost  with  Weight 


Cost  = 1552546.  82  + 78.52  \V 

TV  = 125.  11  x 1013 

R2  = 0. 9624 

EV  = 120.41  x 1013 

UV  = 4.  70  x 10 13 

FQ  = 7936.82 

Cost  with  T -bar 

Cost  = 1827310.07  + 2592502.69  T 

TV  = 125.  11  x 1013  R2  = 0.9117 

EV  = 114.06  x 1013 

UV  = 11.05  x 1013  Fq  = 3199.77 


Cost  with  Structural  Types 

Variable  Name  B^  FQ 


LI  Integral  Blade  593523.  65  4.  17 

L3  Integral  Zee  -82435.  24  0.  08 

Constant  4524667.68 

Note:  L5  had  insufficient  tolerance-level  to  be  included  in  the 

regression  with  an  FQ  = 0.005. 

TV  = 125.  1 1 x 1013  R2  = 0.0166 

EV  = 2.  08  x 1013 

UV  = 123.  03  x 1013  Fn  = 2.  61 
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Table  IV 

Multiple  Regression,  Cost  with  All  Independent 
Variables  Less  T-bar 


Variable 

Name 

3i 

Fo 

W 

W eight 

84.  18 

602. 51 

R 

Rib  Spacing 

-17635. 20 

526. 21 

LI 

Integral  Blade 

-263038. 85 

25.49 

U 

Load 

-0.  14 

2.  60 

L5 

Separate  Jay 

-29905. 29 

0.61 

L3 

Integral  Zee 

-4385. 13 

0.  01 

Constant 

2522902. 36 

Cost 

= 2522902 

.36  + 84 

. 18  W 

- 17635. 20  R - 263038.  85  LI  - 0.  14  U 

- 29905, 

. 29  L5  - 

4385. 

13  L3 

TV  = 

125.  11  x 

1013 

EV  = 

123. 53  x 

1013 

uv  = 

1.  58  x 

10 13 

R2  = 

0.  9874 

TT  - 
x O 

3974. 49 

j-  explanatory  power.  As  rib  spacing  (R)  and  integral  blade  (LI) 

entered,  the  significance  of  U increased  as  evidenced  by  the  increas- 
ing FQ  for  the  coefficient  of  U.  However,  the  computed  F -value  for 

l 

U dropped  from  3250.  54  to  0.81  when  T-bar  entered  the  model.  This 
drop  in  significance  was  attributed  to  the  high  correlation  (value  of  r) 
between  U and  T-bar.  This  relationship  was  further  examined  using 
partial  correlation  techniques. 


Table  V 


Multiple  Regression, 


Cose  with  -All  Independent 


Variables 


W eight 


Variable 

Name 

F0 

U 

Load 

-0.  30 

0.  72 

R 

Rib  Spacing 

-Z3430. 47 

337. 95 

LI 

Integral  Blade 

580850. 26 

86.  94 

T 

T - bar 

3009457. 51 

38.  80 

L3 

Integral  Zee 

-70973. 52 

1.  30 

L5 

Separate  Jay 

-14063. 94 

0.  05 

Cons  tant 

2872221. 05 

Cost  = 287221 . 

05  - 0.  3 U - 23430.  47 

R + 580840. 26 

LI 

-3009457 

. 51  T - 70973.  52  L3  - 

14063. 94  L5 

TV  = 125.  11  x 

1013 

R2  = 0. 9667 

EV  = 120.  94  x 

1013 

UV  = 4.  17  x 

1013 

F0  = 1474. 11 

Step 

Variable 

-ro 

R2 

1 

U 

3250. 54 

0. 9129 

2 

U 

5204.  fc>4 

0. 9458 

R 

187. 36 

3 

U 

7429. 75 

0. 9622 

R 

267.46 

LI 

133.  11 

4 

U 

0.  81 

0. 9665 

R 

339. 29 

LI 

144. 0? 

T 

39.  89 

4 


0. 9665 


PARTIAL  CORRELATION 


To  unmask  any  hidden  relationships  between  or  among  inde- 
pendent variables,  the  partial  correlation  package  in  the  SPSS  pro- 
gram was  used.  Through  the  use  of  partial  correlation,  it  was 
possible  to  examine  the  effects  of  multi-collinearity.  Zero,  first, 
and  second  order  partial  correlation  coefficients  were  obtained  for 
cost  and  each  of  the  independent  variables.  In  each  run,  weight  and 
T-bar  were  the  variables  that  had  the  most  significant  effect  on  each 
of  the  other  independent  variables.  When  weight  and  T-bar  were  con- 
trolled, some  coefficients  which  had  previously  been  negative  became 
positive  and  some  positives  became  negative  (see  Appendix  D).  This 
effect  was  also  noted  in  the  relationship  between  T-bar  and  weight 
where  each  in  the  presence  of  the  other  had  a significant  change  in 
partial  correlation.  For  example,  T-bar  dropped  from  0.95  to  -0.02 
when  weight  was  controlled  and  weight  dropped  from  0.98  to  0.76 
when  T-bar  was  controlled.  The  large  drop  in  the  partial  correlation 
coefficient  of  T-bar  accounts  for  the  seemingly  paradoxical  relation- 
ship obtained  when  cost  was  regressed  with  all  independent  variables 
as  opposed  to  cost  regressed  with  T-bar  alone. 

Partial  correlation  also  explained  the  relationship  of  T-bar 
and  load  (U)  which  was  uncovered  when  cost  was  regressed  with  all 
variables  except  weight.  Recall  that  U showed  a high  degree  of 
significance  until  T-bar  entered  the  equation.  The  zero  order  partial 


for  U with  T-bar  was  0.9981  which  explained  the  effect  T-bar  had  on 
U.  The  zero  order  partial  for  cost  with  U was  0.  9555;  however,  when 
the  first  order  partial  was  taken,  the  coefficient  for  cost  with  U 
became  0.  1349  when  in  the  presence  of  T-bar.  This  was  a significant 
drop  in  partial  correlation  and  accounts  for  the  extreme  drop  in  the 
F-value  for  U when  T-bar  entered  the  regression  model.  Appendix  D 
contains  a graphic  presentation  of  selected  partial  correlation  coef- 
ficients . 

' GRAPHIC  ANALYSIS 

In  order  to  further  analyze  the  relationships  of  cost  with 
weight  and  T-bar,  it  was  necessary  to  plot  the  cost  to  T-bar  and  cost 
to  weight  relationships.  These  graphs  were  generated  on  CREATE 
using  the  CALCOMP  Plotter  (see  Appendix  E for  selected  graphs). 

All  plots  were  obtained  for  a given  rib  spacing  and  structural  type 
over  the  complete  range  of  load  factors.  This  provided  twenty-four 
individual  graphs  for  weight  and  twenty-four  for  T-bar.  To  make  the 
plots  more  meaningful  and  easier  to  interpret,  the  plots  were  drawn 
with  100  points  instead  of  the  observed  thirteen  data  points.  Hence, 
there  was  some  computer  interpolation  between  the  observed  data 
points . 


The  graphs  of  cost  and  T-bar  depict  three  significant  relation- 
ships. First,  the  lowest  cost  occurs  in  a T-bar  range  of  0.25  to  0.50 
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inches.  Second,  T-bar  and  cost  exhibit  a linear  relationship  when 


T-bar  ranges  from  0.5  to  3.0  inches.  For  T-bar  readings  less  than 
0,  5 inches,  there  is  a non-linear  (parabolic)  relationship.  The  third 
relationship  deals  with  the  structural  variables.  The  separate  jay 
consistently  exhibits  the  lowest  cost  and  integral  blade  the  greatest 
cost.  Integral  tee  and  integral  zee  have  the  same  values  and  are,  in 
terms  of  cost,  greater  than  the  separate  jay  and  less  than  integral 
blade.  This  third  relationship  holds  true  for  all  rib  spacings  (26  to 
76  inches).  It  also  holds  for  all  T-bar  values  except  in  the  0.  25  to 
0. 50  range. 

Selected  graphs  of  cost  with  weight  are  included  in  Appendix 
E for  comparison  with  the  cost  with  T-bar  graphs  since  previous 
research  stressed  the  significance  of  the  cost-weight  relationship 
(9:2-10). 
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CHAPTER  IV 


CONCLUSIONS  AND  RECOMMENDATIONS 
INTRODUCTION 

This  chapter  presents  the  conclusions  and  recommendations 
resulting  from  the  research  effort.  The  conclusions  associated  with 
the  research  questions  are  presented  first.  They  are  followed  by 
corollary  conclusions  that  were  developed  in  addition  to  the  central 
research  objective.  A brief  summary  of  the  study  is  presented  and 
then  recommendations  for  future  associated  research  close  the 
chapter. 


CONCLUSIONS 

The  major  thrust  of  this  research  effort  was  to  determine 
the  feasibility  of  developing  algorithms  for  use  with  VDEP  to  optimize 
the  cost  estimation  of  aircraft  structures.  While  the  precise  nature 
of  these  algorithmic  equations  was  not  ascertained,  a significant 
degree  of  insight  was  gained  which  showed  the  feasibility  of  such 
algorithms.  The  research  findings  provide  the  necessary  vehicle  for 
developing  specific  cost  optimization  algorithms.  The  conclusions 
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drawn  from  this  study  are  presented  in  the  same  order  as  the 
research  questions  outlined  in  Chapter  I. 

1.  Based  on  the  statistical  significance  of  the  value 
obtained  through  MLR,  some  relationship  does  exist  between  the 
independent  variables  and  cost.  Variables  found  to  be  the  most 
significant  were  weight,  rib  spacing,  and  integral  blade.  However, 

, 

when  weight  was  deleted,  skin  thickness  (T-bar)  took  over  as  the  *ey 
design  parameter.  Weight  was  not  an  input  design  parameter,  but 


the  result  of  a fallout  from  the  input  parameters  (19).  Hence,  the  net 


effect  is  the  addition  of  two  design  parameters,  T-bar  and  rib  spac- 
ing, to  a suggested  list  of  key  cost  drivers  (9:2-10).  With  these 
design  parameters  the  cost  analysts  will  be  able  to  perform  sensitiv- 
ity analysis  and  make  recommendations  to  the  design  engineers  on 
ways  to  search  for  a least  cost  structure. 

2.  The  relationships  of  load,  T-bar,  rib  spacing,  and  cost 
are  best  depicted  in  the  graphs  contained  in  Appendix  E.  Based  on 
data  collected  for  the  research  effort,  there  exists  a strong  case  for 
linearity.  However,  the  graphs  indicate  that  this  may  not  be  the  best 
interpretation.  The  relationship  shows  signs  of  a much  higher  order 
which  could  be  curvelinear.  Therefore,  a power  function  may  best 
reveal  the  relationship  between  T-bar  and  cost. 

It  is  imperative  to  understand  that  a given  T-bar  and  rib 
spacing  are  two  primary  factors  which  determine  we'ight.  Hence,  the 
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relationshio  between  T-bar  and  rib  spacing  is  an  essential  element 
in  the  development  ox  a successful  cost  optimization  algorithm. 

3.  Based  on  the  statistical  analysis  of  the  data  as  explained 
in  Chapter  III,  the  researchers  believe  that  algorithms  can  be 
developed  which  will,  together  with  VDEP,  assist  design  engineers, 
cost  analysts,  and  program  planners  in  designing  and  selecting  an 
optimum  cost  aircraft  structural  design.  The  principle  result  will 
be  to  strengthen  the  credibility  of  cost  estimations  which  are  vital  to 
the  viability  of  the  overall  design  to  cost  concept. 

COROLLARY  CONCLUSIONS 

The  corollary  conclusions  are  further  results  of  the  basic 
research  effort.  Although  they  were  incidental  to  the  primary  thrust 
of  answering  the  research  questions,  they  add  a great  deal  of  value  to 
the,  original  intent  of  the  study. 

1.  There  exists  some  value  of  T-bar  beyond  which  it  is  not 

feasible  to  search  for  an  optimum  cost.  As  T-bar  reaches  a value  of 

approximately  one  inch,  the  weight  of  the  aircraft  increases  at  an 
> 

increasing  rate.  This  increase  would  require  the  size  of  the  aircraft 
to  be  scaled  upward,  or  alternatively,  a reduction  in  structural  per- 
formance, or  some  combination  thereof. 

2.  Among  the  312  data  values,  minimum  cost  occurred  when 
T-bar  was  in  the  range  of  0.25  to  0.50  inches,  (See  Appendix  E). 
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3.  The  structural  type  selected  was  an  insignificant  factor 
in  determining  cost.  However,  in  the  presence  of  T -bar,  the  integral 
blade  became  an  important  factor  in  the  regression  model  and  was 
significant  in  its  relationship  to  cost.  When  comparing  the  graphs 
for  the  various  structural  types,  it  became  evident  that  the  integral 
blade  consistently  produced  a higher  cost  for  a given  value  of  T-bar. 

4.  As  a general  rule,  when  rib  spacing  increased,  cost 
decreased  while  holding  the  structural  type  constant.  This  relation- 
ship is  depicted  in  the  T-bar  and  cost  graphs  and  leads  to  the  con- 
clusion that  a wing  with  a 56  inch  rib  spacing  can  be  produced  for 
less  cost  than  one  with  a 26  inch  rib  spacing.  In  the  past,  many  air- 
craft structures  have  been  designed  using  a 10  to  20  inch  rib  spacing 
(19).  This  may  have  not  been  the  most  cost  effective  rib  spacing. 
Future  wing  designs  should  consider  rib  spacings  greater  than  20 
inches  to  see  if  cost  can  be  further  optimized. 

SUMMARY 

This  research  effort  attempted  to  uncover  relationships 
between  cost  and  key  design  parameters  in  such  a manner  as  to 
ascertain  whether  algorithms  can  be  developed  which  will  assist 
VDEP  in  designing  an  optimum  cost  aircraft  structure.  By  using  the 
VDEP  program,  it  was  determined  that  some  key  variables  exist 
which  have  a statistically  significant  relationship  to  cost.  Among  the 
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variables  examined,  T-bar  appeared  to  be  the  most  promising  input 
design  parameter  for  use  in  developing  cost  optimization  algorithms. 
Further  study  is  recommended  to  prescribe  the  exact  nature  of  the 
cos t-to-T-bar  relationship  in  order  to  definitize  structural  design 
algorithms . 

RECOMMENDATIONS  FOR  FURTHER  RESEARCH 

1.  Due  to  the  intermixing  of  data  points  in  the  0.  25  to  0.  50 
region  of  T-bar,  additional  analysis  should  be  made  in  an  attempt  to 
spread  out  this  area  and  develop  a better  picture  of  what  happens  to 
cost.  The  area  should  be  enlarged  to  include  smaller  values  of  T-bar, 
even  as  small  as  0.05  inches.  Further  research  in  this  area  will 
assist  in  determining  the  exact  causes  of  the  cost  fluctuations. 

2.  Closely  related  to  the  expanded  study  of  T-bar  is  the 
development  of  a locus  of  points  in  the  lower  ranges  of  T-bar.  This 
locus  of  points  could  be  obtained  by  using  the  cost  with  T-bar  family 
of  graphs.  By  overlaying  the  graphs,  the  locus  of  points  for  lowest 
cost  would  be  plotted.  This  locus  would  then  enable  one  or  more 
equations  to  be  developed  which  would  drive  cost  to  its  minimum 
point. 

3.  The  models  contained  in  this  study  are  highly  simplified 
since  they  are  linear.  An  attempt  must  be  made  to  determine 
whether  higher  order  relationships  exist  among  the  variables  and 
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what  the  higher  order  relationships  are. 

4.  Throughout  this  study,  rib  spacing  (R)  was  the  only 
variable  from  the  group  of  geometric  variables  that  was  taken  into 
account.  Other  variables  affecting  wing  geometry  are  spar  spacing 
and  stringer  spacing.  Therefore  a closer  look  at  the  geometry 
involved  in  designing  an  aircraft  is  in  order.  Spar  spacing  or 
stringer  spacing  may  take  on  additional  significance  and  help  the 
designer  achieve  a more  realistic  cost  optimization  design. 

5.  During  this  study,  data  was  generated  for  the  200th  air- 
craft in  a production  run.  Further  research  should  analyze  data 
generated  by  selecting  the  10th,  40th,  500th,  or  some  other  aircraft 
to  ascertain  the  effect  on  cost. 

6.  Additional  research  should  be  accomplished  to  determine 
which  independent  variable  causes  weight  (W)  to  be  a fallout  rather 
than  a design  parameter.  Partial  correlation  of  all  independent  vari- 
ables with  weight  would  assist  in  uncovering  this  relationship. 
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Algorithm.  "A  set  of  procedures  or  logical  steps  necessary  to  calcu- 
late and  modify  the  input  data  values  to  obtain  the  answers.  It  is  the 
development  of  a logical  approach  to  a problem.  [18:4]." 

Conceptual  phase.  That  phase  of  the  weapon  system  acquisition 
process  which  begins  with  the  determination  of  a needed  operational 
capability. 

This  phase  defines  and  selects  the  system  concepts 
which  warrant  further  development.  The  conceptual 
phase  is  a highly  iterative  process,  with  activities  per- 
formed simultaneously  or  sequentially  as  the  bases  for 
the  acquisition  are  established  by  the  procedural,  fiscal, 
analytical,  experimental,  and  engineering  efforts  accom- 
plished at  DOD  [21:2.  1]. 

Cost  (C).  An  amount  or  quantity  of  money  which  is  given  up  (paid)  to 
purchase  the  specified  goal  or  objective.  All  data  costs  are  given  in 
1971  dollars. 


Design  to  cos t.  A management  concept  wherein  rigorous  cost  goals 
are  established  during  the  early  phases  of  system  development.  The 
control  of  system  costs  (acquisition,  operating,  and  support)  of  these 
goals  is  achieved  by  practical  tradeoffs  between  operational  capability, 
performance,  cost,  and  schedule.  Cost,  as  a key  design  parameter, 
is  addressed  on  a continuing  basis  and  as  an  inherent  part  of  the 
development  and  production  process  (24:2). 
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Design  to  cost  goal. 

A specific  cost  number,  in  constant  dollars,  based 
upon  a specified  production  quantity  and  rate,  estab- 
lished early  during  system  development  as  a manage- 
ment objective  and  design  parameter  for  subsequent 
phases  of  the  acquisition  cycle  [24:2], 


Load  (U'i.  The  amount  of  stress,  measured  in  pounds,  that  is  placed 
on  an  aircraft  structure. 

Multicollinear ity . A condition  that  exists  when  independent  variables 
are  correlated  among  themselves  (14:250).  When  multicollir.earity 
exists,  the  regression  coefficient  of  any  independent  variable  depends 
on  which  other  independent  variables  are  included  in  the  model. 

A regression  coefficient  does  not  reflect  any  inherent 
effect  of  the  particular  independent  variable  but  only  a 
marginal  or  partial  effect,  given  whatever  other  corre- 
lated independent  variables  are  included  in  the  model 
[14:252]. 

Parametric.  A statistical  term  relating  to  an  arbitrary  constant 
which  characterizes,  by  each  of  its  particular  values,  some  particu- 
lar member  of  a system  or  population  (25:1638). 


Performance.  (1)  As  it  applies  to  an  aircraft,  performance  means  a 
prespecified  set  of  parameters:  such  as,  speed,  range,  takeoff  and 
landing  ground  roll,  turn  radius,  etc.  (2)  Pertaining  to  aircraft 
structure,  the  ability  of  a given  structural  type  and/or  material  type 
to  meet  structure  requirements  as  they  contribute  to  overall  aircraft 
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Rib  Spacing  ( R).  Average  distance  in  inches  between  the  stiffening 
element  for  the  box  that  is  the  structural  load  carrying  member 
(19;  15:2-79). 
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Stringer  spacing.  Average  distance  between  the  stiffening  elements 
for  the  skin  cover  of  an  aircraft  wing  (19;  15:2-79). 

T -bar  (T).  Wing  area  spread  over  distance  to  give  an  average  thick- 
ness of  the  aircraft  wing  skin  (19;  15:2-84). 

Weight  (W).  Measurement  in  pounds  of  the  material  usage  to  con- 
struct the  final  product,  a DC- 10  wing. 

Wing.  The  lift  producing  appendage  to  an  aircraft  fuselage.  It  is 
limited  to  the  basic  wing  structure  by  excluding  wing  attachments  such 
as  flaps,  ailerons,  engines,  and  landing  gear. 
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APPENDIX  B 

VEHICLE  DESIGN  EVALUATION  PROGRAM 


1 


i VDEP) 
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The  Vehicle  Design  Evaluation  Program  (VDEP)1  is  intended 
to  be  used  as  a preliminary  design  analysis  tool  to  enable  the  user  to 
rapidly  perform  tradeoff  studies  involving  fatigue,  fracture,  static 
strength,  weight,  and  cost.  The  total  computer  program  is  broken 
down  into  five  modules  (see  Figure  3),  the  program  driver,  vehicle 
sizing,  structural  synthesis,  detail  part  definition,  and  cost  synthesis. 
Output  may  be  selected  from  a complete  and  fully  detailed  version  or 
a summary  version  based  on  the  input  data. 

The  program  driver  calls  the  program  modules  in  proper 
sequence  and  initializes  variables.  It  is  appropriately  titled  in  that 
it  drives  the  total  program. 

The  vehicle  synthesis  module  sizes  the  aircraft,  performs 
a balance  analysis,  distributes  the  area,  and  displays  a planform 
view  along  with  other  pertinent  design  data,  such  as  weight  state- 
ments, center  of  gravity  (CG)  data,  and  general  geometric  data. 

This  program  includes  a curve  plotting  routine  which  allows  the  user 
to  perform  parametric  tradeoff  studies  and  obtain  printouts  for 
further  evaluation.  These  tradeoff  studies  are  possible  at  several 
levels  of  consideration.  For  example,  weight  and  cost  data  can  be 

^For  this  research  effort,  VDEP  was  used  in  the  batch  mode. 
The  interactive  graphics  package  was  not  used.  All  VDEP  processes 
described  herein  were  essentially  the  same. 
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Calls  prograr.  modules  In  the  prope: 
secuer.ce  and  initializes  variables, 


Generates  vehicle  size,  weight, 
and  balance  for  tradeoff  study 
evaluations  during  standard 
preliminary  design  methods. 


Generates  detail  geometry,  loads, 
and  optimum  v;eight  data  for  primary 
structural  components  as  a function 
of  materials  and  construction  type. 


DEFINITIO 


Generates  a list  of  detail  parts  an 
processes  as  a function  of  the 
selected  option.  Computes  actual  an 
material  weight  of  each  detail  part 


Generates  manufacturing  cost  as  a 
COST  function  of  standard  hours  and 

SYNTHESIS  program  cost  and  return-on-invest- 

ment as  a function  of  tooling, 
engineering,  and  manufacturing  cost 


VDE?  Program  Modules 


related  directly  to  key  system  parameters  at  the  vehicle  mission 
level  such  as  payload,  speed,  range,  and  landing  field  length 
requirements.  At  the  vehicle  configuration  level,  data  can  be 


related  directly  to  surface  areas,  spar,  sweep,  taper,  etc.  , and 
fuselage  length,  slenderness,  etc.  At  the  major  component  level, 
comparisons  can  be  made  between  different  materials,  modes  of 
construction,  and  detailed  part  make-up.  Tradeoffs  can  be  made  at 
each  of  these  levels  to  determine  the  overall  vehicle  weight  and  cost 
sensitivities.  This  enables  the  designer  to  refine  the  proposed  air- 
craft design  to  a high  degree  of  detail.  Thus,  engineering  functions 
are  able  to  gain  insight  into  the  cost  effectiveness  of  alternative  air- 
craft system  designs  and  to  determine  the  impact  of  more  detailed 
engineering  alternatives  with  respect  to  any  particular  aspect  of  a 
design. 

The  structural  synthesis  module  provides  detailed  geometry 
loads  and  weight  data  for  the  primary  structural  elements  associated 
with  the  aerodynamic  surfaces  and  the  basic  fuselage  structural  shell. 
Structural  synthesis  provides  a means  of  descriptively  designing 
structural  components  to  fulfill  specified  requirements  of  strength 
and  geometry.  This  program  uses  a multi-station  synthesis 
process.  The  basic  philosophy  of  this  process  is  that  a set  of  struc- 
tural elements  can  be  determined  which  will  satisfy  the  design  cri- 
teria at  each  station  and  that  the  aggregation  of  these  elements  will 
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result  in  a reasonable  representation  of  the  structure.  It  is  the 


structural  synthesis  portion  of  the  total  computer  program  which 
contains  the  subroutine  to  add  or  subtract  material  from  the  struc- 
tural elements  in  an  attempt  to  produce  a minimum  weight  structure. 
The  purpose  of  this  optimization  redesign  process  is  to  change  a 
given  design  to  produce  a lighter  weight  design  while  satisfying  per- 
formance constraints  such  as  minimum  gages,  positive  margins  of 
safety,  etc. 

The  part  definition  program  utilizes  the  output  from  the 
structural  synthesis  module  to  derive  all  the  detail  parts  sufficient 
to  construct  the  complete  assembly  of  the  aircraft.  The  actual  parts 
weights  and  the  weight  of  the  raw  material  to  be  purchased  are  also 
derived  as  part  of  this  module  based  on  the  computed  part  geometry. 

The  part  definition  module  is  coupled  with  the  cost  synthesis 
program  through  the  manufacturing  cost  analysis.  The  manufacturing 
cost  analysis  consists  of  a definition  of  manufacturing  processes 
associated  with  each  part,  the  standard  labor  hours,  and  the  material 
weight.  A list  of  shop  operations  is  called  out  with  each  detail  part, 
and  a series  of  equations  associated  with  each  operation  is  used  to 
compute  the  shop  hours  (standard  hours)  necessary  to  make  the  part. 
By  applying  the  appropriate  labor  rates  to  the  calculated  hours,  the 
direct  and  indirect  manufacturing  labor  costs  are  found.  The  mater- 
ial costs  are  computed  based  on  the  amount  of  material  required  to 
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manufacture  each  part.  Additional  cost  analysis  includes  tooling, 
engineering,  total  program,  and  return-on-investment  costs. 

Tooling  costs  are  computed  as  a function  of  the  number  of 
basic  tool  manufacturing  hours,  initial  and  sustaining  hours  are 
derived  as  a function  of  the  number  of  dissimilar  parts  to  be  pro- 
duced, the  average  number  of  tools  required  per  dissimilar  part, 
and  the  average  number  of  hours  required  to  produce  each  tool. 

Engineering  costs  are  computed  based  on  the  number  of 
manhours  necessary  to  perform  the  various  tasks  associated  with 
the  development  and  production  of  the  aircraft.  Initial  engineering 
hours  are  broken  down  and  distributed  among  the  various  engineer- 
ing disciplines  based  on  studies  made  of  historical  data. 

A learning-curve  approach  is  used  to  derive  costs  of  a 
given  unit  or  lot  as  a function  of  the  first  unit  cost.  Engineering  plus 
the  other  mentioned  costs  comprise  total  program  costs.  Another 
area,  included  for  commercial  operations,  is  a return-on-investment 
cost. 

Each  of  the  five  modules  mentioned  above  has  its  own  input 
data  requirements.  Values  that  are  generated  in  one  module  are  used 
as  inputs  to  other  modules.  This  makes  possible  the  generation  of 
data  within  the  program  which  is  difficult  or  impossible  to  obtain 
during  aircraft  preliminary  design.  Any  time  the  user  desires  to  set 
the  entire  program  back  to  the  original  data,'  it  can  be  accomplished 


■without  having  to  set  each  parameter  to  its  initial  state  individually. 


r 


This  allows  the  user  to  perform  a completely  different  tradeoff 
study  utilizing  initial  input  data  plus  new  changes.* 


The  content  of  Appendix  B is  extracted  from  the  user's  man- 
ual for  VDEP-U,  Computer  Program  to  Assess  Impact  of  Fatigue  and 
Fracture  Criteria  on  Weight  and  Cost  of  Transport  Aircraft,  General 
Dynamics,  Convair  Division,  San  Diego,  CA,  September  1974. 
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APPENDIX  C 

STATISTICAL  ANALYSIS  TECHNIQUES 
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Due  to  the  nature  of  the  data  collected.  Multiple  Linear 


Regression  (MLR)  is  the  appropriate  statistical  technique  to  be 
applied  to  determine  whether  there  is  a relationship  between  cost  and 
the  independent  variables : T-bar,  weight,  structural  type,  and  load 
factor. 

Multiple  Linear  Regression 

The  parameters  of  MLR  are  represented  in  the  population 

equation 

Y = A + BjXj  + B ^ + • • • + BpXp 

where  Y is  the  dependent  variable  (cost),  A is  the  Y intercept,  and  B^ 
are  the  regression  coefficients  of  the  respective  independent  vari- 
ables, X^'s  (15:328).  The  population  equation  is  estimated  from  the 
equation 


y = a + bjX^+b2X2  + . 


• +bpxp  + e' 


where  e refers  to  the  error  between  the  observed  value  of  y and  the 
value  of  y obtained  by  the  regression  equation.  The  objective  of 
regression  analysis  is  to  minimize  cumulative  effects  of  these 
errors  (8). 

The  appropriate  method  for  reducing  these  error  effects  is 
the  method  of  least  squares  (LSM).  LSM  finds  those  values  of  the 
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regression  parameters  which  minimize  the  sum  of  the  squared 


errors,  (e?)  (S’!. 

Several  assumptions  are  required  in  MLR  analysis. 
Assumptions  pertaining  to  the  error  functions  are  referred  to  as  the 
"weak  three"  and  the  "strong  fourth"  (8).  They  are  as  follows: 

1.  e , e are  statistically  independent. 

Xi  J 

2.  Expected  value  of  the  error,  ^E(e„  ))  = 0 when  measured 

xi 

with  respect  to  the  population  regression  line  (A  + BjX^  + B^X?  ~ . . . 

+ BPXP)* 

3.  Variance  of  the  error,  (V(e„  ))  = (variance  of  the 

xi 

population).  This  assumption  is  referred  to  as  the  "homos edas ticity 
assumption"  which  means  'like  variance.  ' 

4.  The  "strong  fourth"  says  that  e is  distributed  normally 

i 

with  a mean  (average)  of  0 and  a variance  of  cr^,  e.  g.  e„  — N(0,cr^). 

xi 

This  assumption  is  necessary  in  order  to  make  statistical  inferences 
concerning  coefficients  a and  b^. 

The  following  three  assumptions  and  limitations  are  also 
necessary  for  the  use  of  MLR: 

1.  The  number  of  sample  observations  must  be  greater  than 
the  number  of  population  parameters  estimated,  n >p. 

2.  The  sample  observations  of  the  independent  variables  are 
linearly  independent. 

3.  Independent  variables  are  observed  or  measured  without 
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error.  Error  is  assumed  to  exist  only  when  measuring  the  dependent 
variable  Y ;8). 

Once  these  assumptions  have  been  made,  the  variables  are 
classified  and  the  analysis  may  proceed. 

Variable  Classification 

The  independent  variables  for  this  research  effort  were  a 
combination  of  two  levels  of  data.  The  information  for  rib  spacing, 
load  factor,  and  weight  was  ratio  level  data.  However,  structural 
types  were  categorical  variables  at  the  nominal  data  level.  For 
these  categorical  variables,  it  was  necessary  to  use  the  method  of 
differences  technique  to  encode  these  variables  for  multiple  linear 
regression  (8). 

The  method  of  differences  designates  one  of  the  levels  of  the 
nominal  variable  as  the  base  level.  The  base  level  should  be  the 
logical  choice  as  the  level  against  which  direct  relative  difference 
comparisons  can  be  made  by  examining  the  variable  coefficients.  In 
the  absence  of  a clearly  defined  advantage  in  selecting  the  variable 
level  which  is  to  serve  as  the  base  level,  a general  rule  of  thumb  is 
to  select  the  level  having  the  fewest  occurrences  among  sample 
observations  (8). 

AFFDL  Structures  Division  did  not  indicate  a preference  for 
the  structural  type  to  be  encoded  as  the  base  level  variable.  There- 
fore, integral  tee  was  randomly  designated  as  the  base  level  of 
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structure.  With  this  base  level  selected,  the  MLR  model  data  was 
recorded  using  a one  (1)  each  time  a given  categorical  variable 
aopeared  in  an  observation  and  a zero  ^0)  otherwise.  P or  example. 


Rib  Soacing 

Int  Blade 

Int  Zee 

Separate  Jay 

Load  F actor 

26. 

1 

0 

0 

374853. 

26. 

0 

1 

0 

1874853. 

26. 

0 

0 

1 

1274853. 

26. 

0 

0 

0 

3674853. 

It  should  be  noted  that  in  the  last  line  of  the  table  above,  the  designated 
base  level  nominal  variable  manifests  itself  as  zeros  in  the  other  nomi- 
nal variables  in  the  data  file.  This  is  a weakness  in  the  method  of 
differences.  "It  is  possible  to  extract  from  the  regression  coefficients 
the  value  of  the  contribution  of  y resulting  from  the  nominal  variable's 
base  level  as  opposed  to  having  no  nominal  variable  information  [8]. 

Analysis 

The  primary  means  for  analyzing  the  data  was  the  Honeywell  63  5 
computer  (CREATE)  at  Wright-Patterson  AF3,  Ohio.  The  data  was 
analyzed  using  the  Statistical  Package  for  the  Social  Sciences  (SPSS), 
version  6.0.  "SPSS  is  an  integral  system  of  computer  programs 
designed  for  the  analysis  of  social  science  data[l5:l]."  SPSS 
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includes  a multiple  regression  capability  in  addition  to  descriptive 


statistics  (15:1). 

SPSS  has  a routine  which  controls  the  order  of  data  entry 
into  the  regression  equation  for  each  of  the  independent  variables 
(15:451).  The  forward  inclusion  scheme  is  the  method  used  to  enter 
the  data.  Independent  variables  enter  into  the  equation  in  their  order 
of  decreasing  net  or  marginal  contribution  to  the  explained  variation 
(8).  A given  variable  is  eligible  for  selection  only  if  its  partial  F 
ratio  is  sufficiently  large,  based  on  criteria  already  established 
within  SPSS.  Partial  F means  the  likelihood  ratio  of  equality  on  the 
test  variable  over  all  groups,  given  the  distribution  produced  by  the 
variables  already  entered.  This  F ratio  is  called  the  F-to-enter 
ratio  (15:453).  The  forward  inclusion  scheme  thus  insures  that  the 
most  significant  variables  enter  the  regression  equation  first. 

In  addition  to  the  regression  equation,  SPSS  output  includes 
the  values  of  three  types  of  variation:  explained,  unexplained,  and 

total.  Explained  variation  (EV)  refers  to  the  component  of  variation 
due  to  or  accounted  for  by  the  regression  (y-"y) ^ (8).  EV  is  also 
termed  the  sum  of  squares  regression  (SSR)  (15:383;  3:465).  Unex- 
plained variation  (UV)  refers  to  the  component  of  variation  due  to, 
unaccounted  for,  or  unexplained  by  regression.  It  is  the  component 
of  variation  about  regression  (8).  UV  is  also  known  as  the  sum  of 
squares  residual  or  sum  of  squares  error  (SSE)  (15:383;  3:4o5). 
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Another  title  for  TV 


Total  variation  (TV)  is  the  sum  of  EV  and  UV. 
is  the  sum  of  squares  total  (SST)  ( 15:383;  3:465).  Thus,  SST  = SSR 
* SSE.  This  concept  is  important  to  help  understand  the  relationship 
that  exists  between  the  dependent  variable  and  the  independent  vari- 
ables . 

Another  output  from  SPSS  MLR  is  the  coefficient  of  determ- 
ination ( R “ ) (15:364).  R“  is  a measure  of  the  extent  to  which  a linear 
multiple  regression  equation  fits  the  observed  data  3:620).  A maxi- 
mum likelihood  estimator  of  the  universe  coefficient  of  determination 

can  be  found  by  dividing  SSR  by  SST  (3:620).  An  explanation  of  the 

2 

test  for  significance  of  R“  is  found  under  the  section  Hyoothes is 
T esting. 

The  regression  coefficients  (B^’sl  are  the  third  major 
section  of  the  SPSS  output.  The  is  important  because  it  represents 
the  expected  change  in  f with  a change  in  one  unit  of  when  all  other 
X^'s  are  held  constant  or  otherwise  controlled  (15:330).  This  makes 
the  regression  coefficients  extremely  important  since  they  determine 
the  key  cost  drivers  among  all  the  independent  variables.  As  with 
R^,  B^'s  have  statistical  hypothesis  tests  for  determining  the  signifi- 
cance of  each  coefficient. 

Hypothesis  Testing 

The  first  test  on  R^  involves  the  overall  regression  equa- 
tion. The  objective  is  to  test  the  statistical  significance  of  the  gross 
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contribution  of  all  the  linear  regressors  simultaneously.  This  test 
indicates  whether  or  not  the  sample  data  evidence  a linear  relation- 
ship between  Y (cost!  and  X , , X?,  . . . , XQ  (8).  The  appropriate 
hypothesis  test  involving  R“  is 
H0:R2  = 0 
H^R2  t 0. 

When  H0  is  true,  Y = BQ  + e.  Hence  one  may  as  well  use  Y to  estimate 
Y.  When  H1  is  true,  Y = (at  least  one  BjX^+e  (8).  The  appropriate 

ccr  ccr 

test  statistic  is  found  by  the  formula  F = -r — — where  p = the 

° ip-1)  (n-p) 

number  of  parameters  in  the  regression  equation  and  n s the  number 
of  sample  observations.  The  test  statistic  iF0)  follows  an  F distribu- 
tion with  p-1  upper  and  n-p  lower  degrees  of  freedom  under  HQ  (7). 

The  F test  is  conducted  as  a one-sided  test.  The  smaller  critical  a 
becomes,  the  more  significant  the  hypothesis  test.  The  more  signifi- 
cant the  test,  the  smaller  the  probability  of  rejecting  a true  hypothesis. 
If  F0  is  greater  than  critical  F (F^),  H0  is  rejected  and  the  overall 
linear  regression  model  is  adjudged  significant.  This  means  at  least 
one  regression  coefficient  is  significant  (8!. 

When  it  is  determined  that  at  least  one  coefficient  is 
significant,  it  is  necessary  to  statistically  test  the  individual  BTs 
to  determine  individual  significance.  The  hypothesis  test  on  the  net 
or  marginal  contribution  of  X^  is 


A 
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H0:Bi  = 0 


H1:Bi  t 0. 


The  test  statistic  is 


it  - SSR  (due  to  X)  x (n-D) 

° sse l8)- 

The  test  on  statistical  significance  of  may  be  conducted 
in  one  of  trwo  ways,  (1)  in  isolation,  or  (2)  with  K other  independent 
variables.  The  conservative  approach  is  taken  by  testing  the  B^ 
values  simultaneously.  This  necessitates  the  use  of  an  equivalent- or 
level  with  which  to  conduct  the  test.  Bonferroni's  approach  computes 
an  equivalent  a = a + K + 1 (where  K = other  independent  variables). 
With  an  a-ievel  of  0.  10,  the  equivalent  ct  is  dependent  upon  the  number 
of  independent  variables  included  in  the  regression  equation  ^8).  If 
F0  > Fc  with  one  upper  and  n-p  lower  degrees  of  freedom  at  an  equiva- 
lent u-level,  the  respective  HQ  is  rejected  and  the  conclusion  reached 
that  the  applicable  B^  is  significant  to  the  regression  equation. 

There  are  three  possible  interpretations  for  a lack  of  statisti- 
cal significance  of  one  or  more  X^'s: 

1.  The  sample  size  is  too  small  to  give  insight  into  the  true 
relationship  among  variables. 

2.  Multicollinearity  is  present.  A solution  is  to  use  a well- 
designed  experiment  where  complete  control  is  exercised  over  all 
independent  variables. 
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3.  does  not  contribute  to  the  regression  equation.  A 
solution  is  to  look  to  a model  without  X If  acceptable,  throw  away 
i8). 


Partial  Correlation  Analysis 

The  basic  formula  for  the  computation  of  partial  coefficients 


■ij  ~ (rik)  irjk> 


where  k is  the  control  variable  and  i and  j are  the  independent  and 
dependent  variables,  r is  the  partial  correlation  coefficient. 

The  extension  of  this  formula  to  more  than  one  control  vari- 
able (that  is,  n + 1 ) is  made  by  replacing  the  simple  correlation  coef- 
ficients (or  zero-order  partials)  on  the  right  side  of  the  equation  with 
the  n^-order  partial  coefficients.  In  this  way  the  preceding  formula 
can  be  used  to  recursively  define  and  compute  each  higher-order 
partial  from  the  previous  one  (15:302-3). 

Partial  correlation  analysis  is  thus  used  as  a tool  to  aid  in 
understanding  and  clarifying  relationships  between  cost,  weight, 
T-bar,  rib  spacing,  load,  integral  blade,  integral  zee,  and  separate 
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SELECTED  GRAPHS 
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This  appendix  is  broken  into  two  sections.  The  first  con- 
tains the  twenty-four  graphs  that  resulted  when  T-bar  was  plotted 
against  cost  for  a given  rib  spacing  (R)  and  structural  type  (Li)  over 
the  range  of  load  (U).  The  second  section  contains  ten  selected  plots 
of  weight  versus  cost  following  the  same  criteria  as  for  T-bar. 
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SECTION  II 
COST  VS  WEIGHT 


R = 7S  .L  = 5 . 1 = 5 
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APPENDIX  F 
COMPUTER  PROGRAMS 


100 


'I 

Multiple  Linear  Regression 


10fr‘jfS,R(SL)  : , 8 , 1 6 ; ; , 1 6 

20$:  IDENT:WP1191, CHRISTENSEN  & EVANS,  7 7A 
30$ : SELECT : SPSS/ SPSS 
4 OR  UN  NAME  ;THES  IS  PROJECT 

50 VAR I ABLE  LIST;R,L5,L1,L3,I5,I3,U,T,C,W 
6 0 VAR  LABELS ; R , RIB  SPACING/ 

70 ;L5, SEPARATE  JAY/ 

80;L1, INTEGRAL  3 LADE / 

90  ;L3, INTEGRAL  ZEE / 

100; 15, ALUMINUM/ 

110;I3, TITANIUM/ 

1 20 ; U , LOAD / 

1 3 0 ; T , T -B AR/ 

140;C,C0ST/ 

150 ;W, WEIGHT 

1 60INPUT  FORMAT ;FREEFIELD 
1 7 OINPUT  MEDIUM; CARD 
180N  OF  CASES ; 3 1 2 

1 90REGRES  SION;  VARIABLES  =*R,L5,L1,L  3,  U,T,  C,W/ 

200;REGR£SSION-C  WITH  R , L 5 , L 1 , L 3 , U , T , W ( 1 ) RESID=0 

2 1 OS  TAT IS  T IC  S ; 2,6,7 

220READ  INPUT  DATA 

230$: S ELECTA: 7 7 A5 7 / THEDATA 

240REGRESSION ; VARIA3L E3 -C , W/ 

250 ;REGRESSION=C  WITH  W(2)  RESID*0, 

260STATI5TICS ; 6, 7 
2 7 OREGRESS ION ; VARIABLES =C ,T/ 

280 ;REGRESSION=C  WITH  T(2)  RESID=0 
290STATISTICS ; 6, 7 

300RECRESSION ; VARI AB LES =C , L 5 , L 1 , L 3 / 
310;REGRE3SION=C  WITH  L5,L1,L3(1)  RESID=0 
320STATISTICS ; 6, 7 
3 3 OF  IN  IS H 
340$ : END JOB 
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Significant  MLR  and  Partial  Correlation 


1 0 V S , ?.  ( 5 L ) : , 3 , 1 6 ; ; , 1 6 

20$:  IDENT:W?U91,  CHRISTENSEN  & EVAN’S,  77A 
30$: SELEET: SPSS/SPSS 

4 OR  UN  NAME  ;THE3 IS  PROJECT  - REGR.  & CORR. 

5 0VARIA3LE  LIST;R,L5,L1,L3,I5,I3,U,T,C,W 
6 0 V A R LA3ELS ;R, RIB  SPACING/ 

7 0 ; L 5 , SEPARATE  JAY/ 

80 ;L 1 , INTEGRAL  BLADE/ 

90;L3, INTEGRAL  ZEE/ 

100 ; 15 , ALUMINUM/ 

110; I 3, TITANIUM/ 

120; u .LOAD/ 

1 3 0 ; T , T -BAR/ 

1 40 ; C .COST/ 

150; W, WEIGHT 

1 60INPUT  FORMAT ;FREEFIELD 
1 70INPUT  MEDIUM; CARD 
ISON  OF  CASES; 31 2 

190REGRZSSION ; VARIABLES =R,L5 , LI, L 3, U, C, W/ 

200  ;REGRESSION  = C WITH  R , L 5 , L 1 , L 3 , U , W ( 1 ) RESID  = 0 

2 1 OS  TATIS  TIC  S ; 2 , 6 , 7 

220READ  INPUT  DATA 

230S: SELECT A: 7 7A5 7 /THEDATA 

2 40RZGR  E S S ION ; VARIA3LE S =R , L 5 , L 1 , L 3 , U , T , C / 

250  ;REGRESSION  = C WITH  R , L 5 , L 1 , L 3 , U , T ( 1 ) RESID  = 0 
260STATISTICS ;2, 6, 7 

2 7 OR EG RE  5 SION; VARIABLES =R , L 1 , C , W/ 
280;RLGRESSION=C  WITH  R,L1,W(1)  RESID=0 
2 9 OS TAT  IS TICS ; 2,6,7 

300PARTIAL  CORR;C  WITH  R , L 5 , L 1 , L 3 , U , T , W BY 

3 1 0 ; R , L 5 , L 1 , L 3 , U , T , W ( 1 , 2 ) 

3 2 OS  TAT  IS  T IC  S ; 1 
33QFINI3H 

340$ : END  JO  3 
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Plots  of  Cost  Versus  T-bar 


10#!?NORM,R(SL) 

20$:  I D E N T : W P 1 191  , AF I T / 5LG  THESIS  PLOTS  - 77A 
30$:M5G2:  1 , SEND  PLOT  TAPE  TO  PLOTTER  (G 0 1 2B , VP  1 1 9 1 ) 

4 0 $ : OPTION:  FORTRAN  , MQMAP 
50$: FORTY : NF  ORM , NLNO 

60  DIMENSION  DUM ( 7 ) , COUNT (13) , X ( 1 3 ) , RX (100) , RY (100) 

70  CHARACTER  LABEL  *7  2 
80  CALL  USTART 

90  CALL  UDIMEN (1 2. , 34. , "CHRISTENSEN  6 EVANS  THESIS\") 

100  10  DO  500  MM= 1 , 3 
110  DO  5 J=1 , 13 

120  READ (16, 1010, END « 999) DUM, COUNT (J) , X(J) 

130  X(J)=X(J)/1000. 

140  1010  FORMAT (V) 

150  5 CONTINUE 
160  YST= (MM-1 ) *1 1 . 

170  YEND-YST+11. 

130  CALL  UPSET ("PRECISION" , 6.  ) 

190  CALL  US  ET ( "SMALL" ) 

200  CALL  UDAREA(0. , 3. 5, YST, YEND) 

210  CALL  UOUTLN 

2 20  CALL  UDAREA (1.5,7.  , YS T+ 1 . , YEND- 1 . 5 ) 

2 30  CALL  US  ET ( "XBOTH" ) 

240  CALL  US  ET ( "Y3  0TH" ) 

250  CALL  UPSET ("XLABEL" , "T-BAR  IN  INCHES\") 

260  CALL  UPSET ("YLABEL" , "COST  IN  THOUSANDS  OF  DOLLARS \ " ) 
270  CALL  USET ("NOLINE") 

280  CALL  US  ET ( "OWN SCALE") 

290  CALL  UWINDO (0. , 3.  , 0. , 12000.  ) 

300  CALL  UPSET ("TICX", .5) 

310  CALL  UPSET ("TICY" , 2000. ) 

320  CALL  UPLOT1 (COUNT, X, 13.  ) 

330  CALL  USET ( "N * " ) 

340  CALL  USET ("SOFT") 

350  CALL  UDOIT ("SETSOFT") 

360  CALL  UMOVE (COUNT (1 ) ,X( 1 ) ) 

370  DO  101  J=i , 13 

380  101  CALL  UPEN (COUNT (J ), X(J ) ) 

390  CALL  USET ("HARDVARE" ) 

400  CALL  USPLIN (COUNT, X, 13. ,RX,RY, 100.) 

410  CALL  US  ET ( "L IN E " ) 

420  CALL  UMOVE (RX(i ) ,RY(1 ) ) 

430  DO  100  J=*l  , 100 

440  1 00  CALL  UPE N ( RX ( J ) , RY  ( J ) ) 

450  READ ( l 7 , 902 ) LABEL 
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4 6 0 

902 

FORMAT ( A 7 2 ) 

470 

CALL 

USET ("DEVICE"  ) 

4 o 0 

CALL 

UMOVE (4 . 25 , YEND-1 . 5 ) 

4 90 

CALL 

UD0IT ( "U?0  1 " ) 

500 

CALL 

UDOIT ( "B  S 0 6 " ) 

510 

CALL 

UPRNT 1 (LABEL, "TEXT" ) 

520 

CALL 

USET ( "VIRTUAL" ) 

530 

5 00 

CONTINUE 

540 

CAL  L 

UERASE 

550 

GO  T 

0 10 

560 

999 

CALL  UEHD 

570 

STOP 

580 

END 

5 90$: LIBRARY: A 1, A2, A3, A 4 
600$: EXECUTE 

610$ : LIMITS : 15 , 40K 

620$: PRMFL: A 1,R,R, GRAPHICS. LIB/GCS/CCS3. 0 

6 30  $ : PRMFL : A 2 , R , R , GRAPHIC S . L 13 /GC S / CALC  3 . 0 
640$: PRMFL : A3  , R , R , AF . L IB / C ALL  1 3 

650$: PRMFL: A4 , R , R , GRADLIB /BATCH 
660$: FFILE: 2 7 , FIXLNG / 8 0 , 3UF S IZ / 8 1 

6 70$ : TAPE : 2 7 , X ID , , , , PLOT -T A? E / WR 
6 80$: DATA: 16 

690$: S ELECTA: 7 7 A5 7 /T HEDATA 

7 00$: DATA:  17 

710$: S ELECTA: 7 7A5 7 /THESIS  04 
720$: EMDJOB 
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1 0 if  v N 0 A .1  , ?,  ( 3 L ) 

20$ : IDENT : N? 1 1 9 1 , AFIT/ SLG  THESIS  PLOTS  - STUDENTS  77A 
3 0 $ : M 3 G 2 : 1 , 5 E N 0 PLOT  TAPS  TO  PLOTTER  (GO  1 2B , WP  1 1 9 1 ) 

40$: OPTION: FORTRAN , NOMAP 
50$: FORTY: NFORM.NLNO 

60  DIMENSION  DUM ( 8 ) , X ( L 3 ) .COUNT (13) , RX ( 1 00 ) , RY ( 1 00 ) 

70  CHARACTER  LA3EL*72 
SO  CALL  US  TART 

90  CALL  UDIMEN ( 1 2. , 34. , "CHRISTENSEN  6 EVANS  THESIS\") 

100  10  DO  500  MM= 1 , 3 
110  DO  5 J=1 , 13 

120  READ (16,  1010, END -99 9 ) DUM, X(J) .COUNT  (J  ) 

130  X(J)=X(J)/1000. 

140  1010  FORMAT (V ) 

150  5 CONTINUE 
160  Y S T = (MM- 1 ) *1  1 . 

170  YEND-YST+11. 

130  CALL  UPSET ("PRECISION"  , 6. ) 

190  CALL  US ET ("SMALL" ) 

200  CALL  UD AREA ( 0 . , 3 . 5 , YS T , Y END  ) 

210  CALL  UOUTLN 

220  CALL  UDAREA(  1 . 5, 7.  , YST  + 1 . , YEND-1 . 5 ) 

230  CALL  US  ET ( "XBOTH" ) 

240  CALL  USET ("YBOTH") 

250  CALL  UPSET  ( "XL  A 3 EL"  , "WEIGHT  IN  POUNDS',") 

260  CALL  UPSET  ("YLA3EL"  , "COST  IN  THOUSANDS  OF  DOLLARS'/’) 
270  CALL  USET ("NOLINE") 

280  CALL  UPLOT1 (COUNT, X, 13. ) 

290  CALL  USET ( "N  * " ) 

300  CALL  USET ("SOFT") 

310  CALL  UDOIT ("SETSOFT") 

320  CALL  UMOVE (COUNT  (1 ) ,X( 1 ) ) 

330  DO  10  1 J = 1 , 13 

340  1 0 1 CALL  U ? E N (COUNT ( J ) , X ( J ) ) 

350  CALL  USET ( "HARDWARE" ) 

360  CALL  USPLIN (COUNT,  X,  13.  ,RX,RY,  100. ) 

3 70  CALL  USET  ("LINE") 

380  CALL  UMOVE (RX(i ) , RY ( 1 ) ) 

390  DO  100  J=1 , 100 

400  100  CALL  UPEN (RX(J ) , RY (J ) ) 

410  READ (1  7, 902  )L A3  EL 
420  902  F0R.MAT  (A  7 2 ) 

430  CALL  USET ( "DEVICE" ) 

440  CALL  UMOVE (4 . 25 , TEND-  1 . 5 ) 

<*50  C.-.LL  UDOIT  ("U?0  1") 
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460  CALL  UDOIT ("BS06") 

470  CALL  UPRNTi (LA3EL, "TEXT") 

430  CALL  USET ("VIRTUAL") 

490  500  CONTINUE 
500  CALL  UERASE 
510  GO  TO  10 
520  999  CALL  UEND 
530  STOP 
540  END 

5 50$  : LIBRARY: A1 , A2, A3, A 4 
560$ : EXECUTE 
5 70$: LIMITS : 15, 40K. 

580$: PRMFL:A1,R,R, GRAPHICS. LIB /GCS/GCS3. 0 
5 90$ : PRMFL : A 2 , R , R , GRAPHIC S . L IB /CC S / CALC 3 . 0 
600$: PRMFL : A3, R, R, AF. LIB /CALLIB 
6 1 0$ : PRMFL : A4 , R , R , GRADLIB / BATCH 
620$:FFILE:27, FIXLNG/80, BUFSIZ/8 1 
630$ : TAPE : 27 , X1D , , , , PLOT -TAPE /UR 
640$ : DATA : 16 

650$ : S ELECTA: 7 7A5 7 /THEDATA 
660$: DATA : 17 

670$: SELECT A: 77A57/THESIS04 
680$ : END  JO 3 
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